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Abstract 

We have recently solved the inverse spectral problem for integrable 
PDEs in arbitrary dimensions arising as commutation of multidimen- 
■ sional vector fields depending on a spectral parameter A. The associ- 

ated inverse problem, in particular, can be formulated as a nonlinear 
Riemann Hilbert (NRH) problem on a given contour of the complex 
A plane. The most distinguished examples of integrable PDEs of this 
type, like the dispersionless Kadomtsev-Petviashivili (dKP), the heav- 
enly and the 2 dimensional dispersionless Toda equations, are real 
PDEs associated with Hamiltonian vector fields. The corresponding 
NRH data satisfy suitable reality and symplectic constraints. In this 
^ , paper, generalizing the examples of solvable NRH problems illustrated 

^ ; in [121 3U H7] , we present a general procedure to construct solvable 

NRH problems for integrable real PDEs associated with Hamiltonian 
vector fields, allowing one to construct implicit solutions of such PDEs 
parametrized by an arbitrary number of real functions of a single vari- 
able. Then we illustrate this theory on few distinguished examples for 
the dKP and heavenly equations. For the dKP case, we characterize 
a class of similarity solutions, a class of solutions constant on their 
parabolic wave front and breaking simultaneously on it, and a class 
of localized solutions breaking in a point of the (x, y) plane. For the 
heavenly equation, we characterize two classes of symmetry reduc- 
tions. 
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1 Introduction 



It was observed long ago [T] that the commutation of multidimensional vector 
fields can generate integrable nonlinear partial differential equations (PDEs) 
in arbitrary dimensions. Some of these equations are dispersionless (or quasi- 
classical) limits of integrable PDEs, having the dispersionless Kadomtsev - 
Petviashvili (dKP) equation [2], [3] as universal prototype example; they arise 
in various problems of Mathematical Physics and are intensively studied in 
the recent literature (see, f.i., [I] - |22]). In particular, an elegant integration 
scheme applicable, in general, to nonlinear PDEs associated with Hamil- 
tonian vector fields, was presented in [5] and a nonlinear d - dressing was 
developed in [H]. Special classes of nontrivial solutions were also derived 
(see, f.i., [13], [16]). 

Distinguished examples of PDEs arising as the commutation conditions 
[I/i(A), L 2 (\)] = of pairs of one parameter families of vector fields, being 
A G C the spectral parameter, are the following. 

1. The vector nonlinear PDE in N + 4 dimensions [25] : 

u tlZ2 - u t2Z1 + (u zi ■ v £ ) u Z2 - (u Z2 ■ V*) U Z1 = 0, (1) 

where U(ti,t 2 , zi, z 2 ,x) G R N , x = (xi,...,x N ) G R N and = (d Xl , ..,d XN ), 
associated with the following pair of (N + 1) dimensional vector fields 

U = d u + \d Zi + U Zi ■ Vs, i = l,2. (2) 

2. Its dimensional reduction, for N = 2 [23] : 

U tx - U zy + (U y ■ V*) U x - (p x ■ Var) U y = 0, 
U G M 2 , x = {x,y), V s =(d x ,d y ), 

obtained renaming the independent variables as follows: ti = z, t 2 = t, x± = 
x, x 2 = y, associated with the two-dimensional vector fields: 

Li = d z + Xd x + U x ■ V^, , v 

L 2 =d t + Xdy + Uy ■ Vaf. l ' 

3. The Hamiltonian reduction ■ U = of ([3]), the celebrated second 
heavenly equation of Plebanski [21] : 

o tx -e zy + e xx e y y-e 2 x =o, e = 9(x, y ,z,t) ei, x, y ,z,teR, (5) 
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describing self-dual vacuum solutions of the Einstein equations, associated 
with the following pair of Hamiltonian two-dimensional vector fields 

L\ = d z + Xd x + 6 xy d x — 6 xx d y , ^ 

L2 = 9 t + Xdy + Oyyd X ~ X ydy . 

4. The following system of two nonlinear PDEs in 2 + 1 dimensions [25] : 



U X t ~\~ Uyy ~\~ {%LU X ^ X ~\~ V X 11 X y VyU XX 0, 

V X t ~\~ Vyy ~\~ UV XX ~\- V x V x y VyV xx 0, 

arising from the commutation of the two-dimensional vector fields 

Li = d y + (A + v x )d x - u x d x , 

L 2 = d t + (A 2 + Xv x + u- v y )d x + (-Xu x + u y )d x , 

and describing a general integrable Einstein- Weyl metric [26J. 
5. The v = reduction of (|7j), the dKP equation 



(7) 



(8) 



(u t + uu x ) x + u yy = 0, u = u(x,y,t) el, x,y,teR, (9) 

(the x-dispersionless limit of the celebrated Kadomtsev-Petviashvili equa- 
tion [27J), associated with the following pair of Hamiltonian two-dimensional 
vector fields [3 E] : 



Li = d y + Xd x - u x d x , 

L 2 = d t + (A 2 + u)d x + {-Xu x + u y )d X: 



(10) 



describing the evolution of small amplitude, nearly one-dimensional waves in 
shallow water [28] near the shore (when the x-dispersion can be neglected), 
as well as unsteady motion in transonic flow [2] and nonlinear acoustics of 
confined beams [5]. 

6. The u = reduction of © EZ]: 

V X t+Vyy =VyV XX ~V X V X y, V = V (X , , t ) G R , I, ]/, t £ I, (ll) 

associated with the non-Hamiltonian one- dimensional vector fields [IS] 

L 1 =d y + (X + v x )d x , . , 

L 2 = 9 i + (A 2 + A^-^)9 a; . 1 j 
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7. The two-dimensional dispersionless Toda (2ddT) equation [221 ED]: 

<Pc^ = {e^) t , = 0(Ci,C 2 ,t) (13) 

(or <^CiC2 = ( eip )tt ' V 9 = 0f)> associated with the pair of Hamiltonian vector 
fields 0: 

L x = d Cl + \e%d t + (-A( C ^) t + %i) Ad A , 

L 2 = % + A-^ft + (A-^e^), - %i) Ad A , (M) 

describing integrable heavens [311 [32] and Einstein - Weyl geometries |33j, 
[31], [35]; whose string equations solutions [10] are relevant in the ideal Hele- 
Shaw problem [3"B]-[4"0"]. 

The Inverse Spectral Transform (1ST) for 1-parameter families of multi- 
dimensional vector fields, developed in [23] , has allowed one to construct the 
formal solution of the Cauchy problem for the nonlinear PDEs ([3]) and (jSJ) 
in [23], for equations flTJ and in [25], for equation f TTTj) in [UJ and for the 
wave form (e^*)t = + 0^ of equation ( fl3l) in [42]. This 1ST, introducing 
interesting novelties with respect to the classical 1ST for soliton equations 
[431 128] , turns out to be, together with its associated nonlinear Riemann - 
Hilbert (NRH) dressing, an efficient tool to study several properties of the 
solution space of the PDE under consideration: i) the characterization of a 
distinguished class of spectral data for which the associated nonlinear RH 
problem is linearized and solved, corresponding to a class of implicit solu- 
tions of the PDE (for the dKP and 2ddT equations respectively in [UJ and 
in [12], and for the Dunajski generalization [16] of the heavenly equation in 
[43] ) and for equations (ITT]) and ([3]) in [1?]; ii) the construction of the long- 
time behaviour of the solutions of the Cauchy problem (for the dKP, 2ddT 
and heavenly equations respectively in [37], [H] and [UJ); iii) the possibil- 
ity to establish whether or not the lack of dispersive terms in the nonlinear 
PDE causes the breaking of localized initial profiles (for the dKP, 2ddT and 
heavenly equations respectively in [UJ, in [12] and in [47]) and, if yes, to 
investigate in a surprisingly explicit way the analytic aspects of such a wave 
breaking, as it was done for the dKP equation in [UJ. Recent results on 
integrable differential constraints on the hierarchy associated with the non- 
linear system (J7J) and their connection to the associated NRH problems can 
be found in [48J. 

In this paper, generalizing the examples of solvable NRH problems illus- 
trated in [1211111117], we present, in §2, a general procedure to construct solv- 
able NRH problems for integrable PDEs associated with Hamiltonian vector 
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fields, allowing one to construct implicit solutions of such PDEs parametrized 
by an arbitrary number of real functions of a single variable. In §3 we illus- 
trate this theory on few distinguished examples for the dKP equation, includ- 
ing the similarity solutions, a class of solutions constant on their parabolic 
wave front and breaking simultaneously on it, and a class of localized solu- 
tions breaking in a point of the (x, y) plane. In §4 we briefly apply the theory 
to the heavenly equation, constructing their similarity solutions. 

Since the theory will be illustrated on two basic examples of PDEs asso- 
ciated with Hamiltonian vector fields: the dKP and the heavenly equations, 
in the remaining part of this introductory section we summarize their NRH 
dressing formalisms. 

1.1 NRH dressing for dKP |25].g44] 

Consider the vector nonlinear RH problem on the real line: 



^+ = 7e 2 (Vr,V 2 ~), Ael, (15) 
or, more shortly, 

#+(A)=£(#-(A)), (16) 

where TZ j(s\, S2), j — 1, 2 is a given pair of complex different iable functions 
of two arguments, and the solutions ^(X) = (V>i (A), ^j(A)) G C 2 are 2- 
dimensional vector functions analytic respectively in the upper and lower 
halves of the complex A plane, with the following asymptotics, for |A| 3> 1 
in their analyticity domains: 



#(A) = -XH -Xy + x- 2tqi 1) + E V- 

n>l 

„(!) _ „(") 



(17) 



A 1 {-< A ' 

n>2 



where, f.i.: 



« = - yu + 2td- l u y , qf ] = d~\yu) y - 2t<9 x 7 2 u 



mi ■ 



q?>=u, q 2 2) = -d^u y , qf = -u 2 /2 + d- 2 u yy , (18) 
= \dx l {u 2 ) y + 2d- 1 (u x d- 1 u y ) - d- 3 u yyy . 
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Then, if the nonlinear RH problem (j!5p is uniquely solvable, together with 
its linearized version 

M ifc (v>r,^ = ^%^, j,k = i,2, (19) 

then ipf 2 are solutions of the following pair of vector field equations: 

Lyb = ip y + (A + v x )4) x - u x ip x = 0, 

L 2 tp =i>t+ (A 2 + Xv x + u - v v )4j x + (-Xu x + Uy)lfj\ = 0, 

where 

(i) 

u = q 2 , 

(1) (1) o+ (2) 

v = - q \' - yq ± - 2tq\ \ 
and u and v solve the system ([7j) of PDES. 



(20) 



(21) 



2. Basic non-differential reductions. There are few basic non-differential 
reductions of the above NRH problem and, correspondingly, of the integrable 
system (j7j). 

Rl. If TZo(s 1,82) = s 2 , then 1P2 = 4>2 = A, the nonlinear RH problem 
becomes scalar and one obtains equation (jTTI) . the u = reduction of the 
system ([7]), associated with the non Hamiltonian Lax pair of one - dimen- 
sional vector fields (j!2p . 

R2. If the transformation (si,s 2 ) —> (7£i(si, s 2 ), TZ 2 (si, s 2 )) is canonical: 

{7^-1, 7^ 2 }( Sl , S2 ) := ^lsi^2 S2 - T^lsz^s! = 1> ( 22 ) 

one obtains the v — reduction of the system (J7J), the celebrated dKP 
equation ([9]), corresponding to the Hamiltonian Lax pair 

ipy + Xip x - u x ip x = ip y + {H 2 , ip}(\,x) = 0, (23) 

i) t + (A 2 + u)i) x + {-\u x + %)V>a = + {#3, ip}(\, x ) = 0. (24) 



for the Hamiltonians 



A 2 A 3 
H 2 = — +u, H 3 = — + \u- d x x u y . (25) 
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We remark that, in this reduction, the eigenfunctions ipf,ipf are canonically 
conjugated: 

{^^f}(x, x) = l (26) 

We also recall that [6j [5] the dKP equation is the first member of the dKP 
hierarchy 

H mtn - H ntm + {H n , H m } {x ,x) = 0, m ^ n, m, n > 2, 
if B = + , n>2, 

corresponding to the Hamiltonian vector field equations 

i> tn + {H n ,i/;} (x>x) = 0, (28) 

where t 2 = y, t 3 = t and {jp)+ is the non negative (principal) part of the 
Laurent expansion of ip at A = oo. In particular, if m = 2 in (|27|) . one 
obtains the sub-hierarchy 

nu tn + (Res^n), = 0- (29) 

where Res(g) is the coefficient of A" 1 in the Laurent expansion of g(X) at 
A = oo. 

R3. If 

£(£(C)) = C, VCGC 2 , (30) 
then one obtains the reality constraint 



u,veR, ^+(A)=^-(A), j = 1,2. (31) 

From the integral equations characterizing the solutions of the NRH prob- 
lem ( Tl5|) . and from the definition u = q£ , one obtains the following spectral 
characterization of the solution u: 

u = F(x - 2ut,y,t) e R, (32) 
where the spectral function F, defined by 

F (£, y,t) = -J ^R 2 (vrr(A; £, y, t), rc^X; £, y, tj) , 

^(Ci,C 2 ) = + ^(Ci,C 2 ), j = 1,2, 
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is connected with the initial data u(x, y, 0) via the direct spectral transform 
developed in [23]. Equation ( 13"2"j) defines the dKP solution implicitly , due to 
the presence of the (x — 2ut) term as argument of F, and describes the wave 
breaking features of localized solutions of dKP. 

Evaluating the integral equations characterizing the above NRH problem 
in the space-time asymptotic region 



x +u=^ % = £ + vit, y = v 2 t, ^ 



£,-2ut,v h v 2 = 0(l), v 2 ^0, £>1, 

2 

one obtains the following longtime behavior 

u = 7t G ( x + V Tt- 2uh ii) +0 {-7t l (35) 

where 

G(^,r/) = -^r J dfxR 2 (^ - /i 2 + ai{ti;Z,r}),-Ti + a 2 (/i;C,^)), (36) 

K 

and ai t2 (fi; £,,rj) are the unique solutions of the nonlinear integral equations 

k V (37) 

-r) + a 2 (ii';£,r))j, j = 1,2, 

giving a description of the wave breaking of small initial data in the longtime 
regime as explicit as for the case of the Riemann Hopf equation u t + uu x = 
0. In particular, one shows that small and localized initial data will break 
asympotically in a point of the (x, y) plane and the breaking details are given 
explicitly in terms of the initial data 



1.2 NRH dressing for heavenly [[231, M\ 
Consider the vector NRH problem on the real line: 

^ + (\) = tf-(\) + R(tf-(X),X), AgM, (38) 

where ip + (X), if}~(\) G C 2 are 2-dimensional vector functions analytic in the 
upper and lower halves of the complex A plane, normalized as follows 

H^)=( X y Z X x Z t )+0(X- 1 ), |A|»1, (39) 
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and the spectral data R((, A) = C2, A), i? 2 (Ci, C2) A)) G C 2 , defined for 

( G C 2 , A G IR, satisfy the following properties: 



TZ(JZ(C, X), X) = C, VC G C 2 , reality constraint, ^g-j 
{Hi, 1^2}^ = 1, symplectic constraint, 

where 1Z((, A) := C + ^(C>-\)- Then, assuming uniqueness of the solution 
of such a NRH problem and of its linearized version, it follows that if} are 
canonically conjugated solutions ({ipf, ^f}^^) — 1) °f the linear problems 
Li^^ = 0, where L 12 are defined in ([6]), and 

^ )=F>,|,,^)6M 2 (41) 

is solution of the heavenly equation ([5]), where 

dX 



F(x, y, t) = / ^—R{^Pi 0; x, y, z, t), -02 (A; x, z, t), a) . (42) 



As a consequence of equations Li^ip^ = 0, it follows that, for |A| 3> 1: 

^± = x - Xz - OyX- 1 + e t x- 2 + o(A- 3 ^ 



tjjf = y - Xt + ^.A- 1 - ^A~ 2 + 0(\- 



3> 



(43) 



2 Solvable vector nonlinear RH problems 

Since basic nonlinear PDEs like the dKP, heavenly and 2ddT equations are 
solved via NRH problems whose data satisfy the symplectic and reality con- 
straints, in this section we present a general procedure to construct solvable 
vector NRH problems satisfying such constraints. 

Consider an authonomous Hamiltonian two-dimensional dynamical sys- 
tem with Hamiltonian 

H(x) = H(S(x)) (44) 

where x = (q,p) are canonically conjugated coordinates, %{•) is an arbitrary 
entire function of a single argument and S(x) is an entire function of the 
coordinates, corresponding to the equations of motion 

f =^)(_° 1 J)v^(x). (45) 
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Introducing action - angle variables in the usual way: 
J = ±$p(q,H)dq, =*► H = H{J), 



9-9 = u(J)H'(S)(t - t ) = H'{S) J dp(q d H / J)) dq', (46) 

go 

,.( j\- 9H(J) 

the solution can be found inverting the quadrature (jlSj) : 

x = V(9-9 ;x ,J) (47) 
where the evolutionary diffeomorphism T> : xo — )■ x is symplectic: 

{V u V 2 } {qo , P0) = l (48) 

and satisfies the properties: 

x = V(9 -9;x,J), 

V(9 2 ; V{9 x] x , J),J) = T>{9 X + 9 2 ; x , J), (49) 
V(0;x o , J) = x . 

Consequently, we have the relation 

V{9 X - 6; x, J) = V(9i - 9 ; x , J) (50) 

for any intermediate angle 9 X , and this relation holds true, at least locally, 
also in the complex domain. 

2.1 Solvable NRH problems 

Identifying the solution x(t) at time to and at time to + 1 with, respectively, 
if-(X) and ^ + (A), 

i) equation ( 14"?]) becomes the two-dimensional vector NRH problem 

J + = v (u(j$-))u\s$-)y, r , j$- )) = n ($-) (si) 

connecting the (— ) and (+) vector functions through the canonical transfor- 
mation (122]); 

ii) since Six) (as well as J(x)) is an invariant of the dynamics: S(xq) = S(x), 
it follows that S(ip~) = S(ip + ) is an invariant of the NRH problem ( 15"TT) . 
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Therefore 5 , (^ ± ) define an entire function, given by their polynomial part 
W(X): 

S(J-(\)) = S(J + (\)) = W(\). (52) 

In addition, since S is a real function of its arguments, the reality constraint 
( 13 ip implies that 

W(X) = W(X). (53) 

iii) Since also T> is a real function of its arguments and, from (I49p . 

V(z;V{-z;(,J),J)=C, (54) 

the reality constraint (|30|) implies that V(—z;(,J) = V(z;(,J); it follows, 
from (??), that 

«'(•) = </(•), (55) 
where / is an arbitrary real function of a single argument, for dKP, and of 
two arguments (the invariant (|52|) and A), for heavenly. 

iv) Defining the upper and lower functions 

9± w s h I y^±m) (™wrm))mm)) m, (5 6) 

K 

so that iuf = 9 + — 6~ , becomes 

P(-^ + ; J(^ + )) = V", )), (57) 

and provides the solution of the NRH problem, at least if T> is an entire 
function of its arguments. Indeed, in this case, the left and right hand sides 
of ( IBTj) are equal to their polynomial part E{X) at A ~ oo: 

V(-6 ± ;j ± J J(j ± )) = E(\) (58) 

and these are algebraic equations for ^(X). 

Remark 1. Solvable NRHs of this type allow one to construct solutions of the 
dispersionless PDE characterized by differential reductions of the integrable 
PDE. To show it, let us consider, for concreteness, the dKP case. Let g(A) be 
a solution of the NRH problem, eigenfunction of the spectral problem (123j) . 
and characterized by the formal Laurent expansion: 

N 

g(X) ~ 9nX n , |A| » 1. (59) 
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Then the coefficients g n of A™ satisfy the recursion relations 



9n x — 0' 9n v + 9n-i x — 0, 

9n y + 9n-i x ~{n+ l)u x g n+1 = 0, n < N. 



(60) 



Since W(\), defined by ( )52|) . is a polynomial (in A) solution of ( )23|) . then 
Res oc W / (A) = W_i = and equation ( 160]) for n = yields: 



In addition, equations ( 160]) for n < imply, as it has to be, that W n = for 
n < — 1. Equation (16T]) is the differential constraint satisfied by the solutions 
of dKP associated with the solvable NRH problem. 

Remark 2. Starting with any In - dimensional Liouville integrable Hamil- 
tonian system, a generalization of the above construction allows one to ob- 
tain a solvable 2n-dimensional NRH problem ip + (X) = lZ(ip~(\)), where 
1Z : ip^ ip + is a symplectic transformation. 

Remark 3. Also integrable symplectic maps from W 1 to W 1 generate, following 
the same procedure, solvable n dimensional vector NRH problems. 

2.2 Increasing the richness of the solution space 

In the previous section we showed how to construct solvable NRH prob- 
lems whose spectral data, satisfying the reality and symplectic constraints, 
have the freedhom of an arbitrary real function of a single variable, which 
parametrizes the corresponding space of implicit solutions of the dKP, heav- 
enly and 2ddT equations. In this section we show how it is possible to increase 
considerably the richness of the space of implicit solutions. For concreteness 
we refer to the dKP equation, but the same considerations apply to all real 
PDEs associated with Hamiltonian vector fields. 

Let ijjf 2 be solutions of the Hamiltonian vector field equations (123]) . (]24]) , 
satisfying the asymptotics ( 117]) . then arbitrary differentiable functions of 
{^iii^f) are a ^ so solutions of (]23]) . (124]) . In the first step of our construc- 
tion we define new functions of (ipi , ip 2 ) as follows: 



Woy - u x W x = -W. lx 



0. 



(61) 



4 



(i)± 



4 



(i)± 



(62) 
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where ai G R and is an arbitrary real function of one variable. The trans- 
formation (I6"2"j) is clearly invertible and, most important for our purposes, is 
symplectic and preserves the reality constraint ( IBTj) . Therefore Vi 2 are also 
canonically conjugated solutions of the vector field equations (I23"]) . fl24")) satis- 
fying the reality constraint: {^| 1 ' )± , 1^2^} = 1) ^j^M — (A), J = 1)2, 
parametrized by the arbitrary real function of a single argument. In the 
second step of this construction, we use the solutions V , i^2 ± ^° construct new 
functions as follows 

Vi? )± = ^ 1)± , 

,(2)± ,ri)± . ,f 2 i ^ /C 1 ^ (^^) 



where a.2 G M and / < - 2 - ) is another arbitrary real function of one variable (notice 
that, while in (|62|) the transformation is identical wrt the second eigenfunc- 
tion, in ( 163]) the transformation is identical wrt the first one). Therefore 

( 2)i 

we have constructed the family tp\ 2 of canonically conjugated solutions 
of the vector field equations f l23|) . fT24|) satisfying the reality constraint, and 
now parametrized by the two arbitrary real functions f^',f^ 2 ' of a single 
argument. This procedure goes on, without limitations, alternating transfor- 
mations of the type (l62|) and (I63I) . which are identical respectively wrt the 
second and first eigenfunctions. At the m th step, one constructs canonically 
conjugated solutions ip^ ± of (!23l . fl24l) satisfying the reality constraint: 



ty} m)± ,^ m)± } = 1, )+ (A) = ^ (mK (A), j — 1,2 (64) 

and parametrized by m arbitrary real functions . . . , of a single ar- 
gument. 

Now we consider any of the solvable NRH problems, for the family of 
eigenfunctions ± ! f° r instance, the NRH problem discussed in §3.2: 

whose normalization easily follows from the definition of i/>i2 ^ ^ n t erms of 
Vt~2> an d from the asymptotics (TPTj) . Then the solution of this exactly solvable 
NRH problem and the corresponding family of exact implicit solutions of dKP 
are parametrized by m + 1 arbitrary real functions /, /W, . . . , /( m ) of a single 
argument. 
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3 Examples for dKP 

In this section we consider some basic examples of solvable NRH problems 
associated with dKP. 

3.1 Example 1. The invariant ipiip^ an( ^ the similarity 
reduction 

If S(q,p) = qp, then equation ( 14"T|) reads 

q(t) = q,e W{s ^- t0 \ p(t) = poe^'^^; (66) 
and corresponds to the NRH problem 

= ^- e W^) ) ^+ = 1 ^- e -#r*2) i Aei, (67) 

satisfying the symplectic and reality constraints. Then the invariance equa- 
tion ( 1B^|) becomes: 

tfj+^+ = = -t\ 3 - y\ 2 + (x- 3ut)X - 2yu + 3t d~ x u y = W(X), (68) 

and the NRH problem linearizes and decouples: 

^+ = ^-e^ w ^\ = ^ 2 "e- i/(w(A)) . (69) 
Equations f )58|) become 

^+ e i(-y /+ (A) = ^,- e i(-y/-(A) = jEj(a); J = 1; 2 , (70) 

where the analytic functions / ± (A), defined by 

fH^^^Ly^fiWiX')), (71) 
exhibit the following asymptotics 

/ ± (A)~^E(A"- 1 /)A-", 

n>l (72) 

(AV> = i/ H A-(W(A))dA, 

if / decays faster than any power. It follows that 

(A) = -tA 2 - (y + t < f >)A + x - 2ut - y < f > - . , 
t(<A/>+i</> 2 ), E 2 (X) = X- < f >, lMJ 
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implying the following explicit solution of the NRH problem: 

rl>f = E j (\)e«-y i+1 f + W, j = 1,2. (74) 

A characterization of the corresponding solutions of dKP is obtained observ- 
ing that, since W in (I6"5j) depends on the unknowns u and d~ x u yi isolating 
the 1/A term of equations (1741) for |A| ^> 1 and comparing them with the 
asymptotics (fT7j) . (fTHj) . we obtain the algebraic system 

qf ] = 2td~ l u y -yu=-{x- 2ut) < f > +y(< f > 2 /2+ < \f > + 
t(< \ 2 f > + < f >< \f > + < f > 3 /6)), (75) 
= u =< \f > — < f > 2 /2. 

for the unknowns u and d~ x u y . The constructed solutions of dKP correspond 
to the following differential reduction (ISTj) : 

3tu t + xu x + 2yu y + 2u = 0. (76) 

Substituting its general solution 

u = r 2/3 A(x',y') , x' = -^, 1/=^, (77) 

into dKP, one obtains the similarity reduction of dKP: 
x' 2 

A x > + —A X ! X , + -y'A x , y t - Ay, yl - ( A A x > ) x , = 0. (78) 

Therefore the algebraic system (175]) characterizes these similarity solutions 
of dKP. 

3.2 Example 2. The invariant i/j^ + cli/j^ 

If S(q,p) = q + ap, where a is a real parameter, then equation (1471) read 

= q + a-H'(5)(t - t ), p(t) =p - W(S)(t - t ), (79) 

becoming the NRH problem 



15 



satisfying the symplectic Q and reality (13"UI) constraints. 
Due to the invariance equation 



if>£ + aip$ = + a ^2 = -(y-a)X + x- 2ut = W{\), (81) 

the NRH problem linearizes and decouples: 

1>t = fo+iaf(W), tfj+ = ^2-^f(W), AgR, (82) 

and equations f[5"8l become 

^+(A) - iaf + (\) = - iaf-(X) = E 1 (X), 
^+(X) + ^f + (X) = ^ + if-(X) = E 2 (X), (83) 
(A) = -tX 2 -yX + x-2ut, E 2 (X) = X. 

Then the solution of the NRH problem reads: 

ifjf = -t\ 2 -yX + x-2ut + iaf* (A) , (84) 

il>t = \-if ± (\). (85) 

Using equations (I8"4"j) . (l8"5]) and (fT7j) . it is possible to express the coefficients 
q^l of the asymptotic expansions in terms of the spectral function f(W) in 
the following way: 

g W = - a < X n ~ l f >, n > 1, (86) 

g W =< A™- 1 / >, n > 1. (87) 
Since W in flSTl) is function of = u only, equation f|8T|) for n = 1 

^M i+ t- 2 4 (88) 

characterizes the family of solutions associated with the above NRH problem, 
where 

F{z) = -L / /(-/i 2 + *)t^. (89) 



2vr 

Equation (1881) describes a family of solutions of dKP, constant on the parabola 

% + J t = f , (90) 

16 



and breaking simultaneously in all points of it. Therefore the above exact 
solutions do not describe the breaking of a localized wave in the (x, y) plane, 
typical of the dKP evolution of localized waves jH] . In the next subsection 
we show how to generalize this solvable case to construct exact solutions 
describing the breaking of a localized wave in a point of the (x, y) plane. 
The corresponding differential constraint (IBTj) reads 

(y — a)u x — 2tu y = 0. (91) 

Indeed, substituting the general solution 

u = <pfot), Z = x+ {V ~ t a) \ (92) 

of ( 191]) into the dKP equation, one obtains the equation (pt+tp/ (2t)+ipip£ = 0; 
at last, setting tp = t _1 / 2 t>(£, r), r = 2\/i, we finally derive the Burgers - 
Hopf equation v T + vv% = 0, whose general solution reads v = F(£ — vr), 
where F is an arbitrary function of a single argument. Going back to the 
original variables, we recover 



3.3 Example 3. The invariant ^ + a(^) n 

If S(q,p) = q + ap n , where a is a real parameter, then equation ( 1471) read 

q(t) = q + ap n -a(p -'H'(S)(t-t )) n , p(t) = p - U'{S){t - t ), (93) 

becoming the NRH problem 



V> 2 + = ^ ~ if&i + aVO, A G R, 

whose data satisfy the symplectic ([9]) and reality (1301) constraints. 
Due to the invariance equation ([321) 



(94) 



ipt + aipt™ = ipi+a^ n = -t\ 2 -y\ + x-2ut + a ((^)") + = W(X), (95) 

the NRH problem linearizes and decouples and equations (1551) become (1551) 
and 

^+(A) + if+(\) = V£ + if- (A) = A, (96) 
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Then the solution of the NRH problem reads: 



1>i = W(\) - a^ n , (97) 
ifit = \-if ± (\). (98) 

Since now W depends on the (n — 1) unknowns u, q% , n = 2, . . . ,n — 1, 
the system of (n — 1) equations (!H7j) . for n = 1, . . . , n — 1, is an algebraic 
system characterizing a family of implicit solutions of dKP parametrized by 
the arbitrary real function / on a single variable. 

The corresponding differential constraint (ISTj) reads, due to fl29|) . 



(^(ipf + aipf n ) ^ = yu x - 2tu y + anu tn = 0, (99) 

where u tn is the n th flow of the dKP sub-hierarchy (1291) . 
For n = 2, the invariant 

^ + aipf 2 = -t 2 \ 2 -y\ + x- 2ut 2 = W{\) (100) 

where r 2 = t — a, is equivalent to f lHTj) . up to trivial shifts of y and t; therefore 
the corresponding solution of dKP is essentially the same as that in §3.2. 

Now we show that, if n > 2 and even, this family of solutions describes 
the evolution of a localized two-dimensional wave breaking at a point of the 
(x,y) plane. To do it, it is convenient to view ( 194|) as a NRH problem on 
the real line in the shifted variable //: /i = A + y/2t, t ^ 0, and rewrite the 
asymptotics (TT7I) in terms of /i: 



n>l 

~(n) 

n>l 



(101) 



(102) 



where, for instance: 

HI — u £ a rii HI — 2t €, W 

and the new space variables £ and r\ are defined by 

2 

£ = £ + — , 77 = —. (103) 

s 4t' 7 2t v ; 
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In addition, equation (157)1 is replaced by 

= (^f) = ^J ^fiWi^W. (104) 

R 

where W(fi) = (^f + GM/>^ n )+ is now the non negative part of the Laurent 
expansion, in the \x variable, of the invariant. 

If n = 3, then 

W(fj) = afi 3 - r 3 fi 2 + 3a(r] 2 + u)fi + £- ar] 3 - 2ur 3 + 3qf\ (105) 

where 

r 3 = t + 3ar]. (106) 

Since W depends on the two unknowns u, q\ , the system of two equations 
( I104p . for n — 1, 2, is an algebraic system characterizing a family of implicit 
solutions of dKP parametrized by the arbitrary real function / of a single 
variable. 

It is remarkable that, in the longtime regime and for \a\ <C 1, it is possible 
to obtain the following explicit asymptotic formulas for such a solution: 



u ~ -$=F{£ — an 3 — 2ut 3 i, 

7/ = 0(1), £-ar? 3 - 2wr 3 = 0(1), t > 1, \a\ < 1, 



where F is defined by (1891) . 

To show it, we first observe that the condition \a\ «C 1 implies: 

J R f(W(v)W ~ / r [/(W(°)(m)) + a/' = 

it/G^M) + a/' (WW^JWSV)]^ ~ ( 108 ) 
J M f(WM(/i) + aWV(ti)W, 

where = W<® + aW^\ with 

^(°)(^) = ^f }± ) = + f - 2itt, 

W^ 1 )^) = (V4 1)±3 ) = /i 3 - 3r^ 2 + 3(r/ 2 + u)/i - - 6r]u + 2,qf\ 

(109) 

and Wev (a*) is the even part of W^ 1 ^//): 

= -377^ - r/ 3 - 6r?M + Zqf \ (110) 
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Second, the conditions f > 1, \rj\ = 0(1) imply that r 3 ^> 1 and suggest 
the change of integration variable // = yffi[i in (11081) . Consequently u = 
0(1/ ^J%) = 0{l/y/t)\ the same change of variables in the integrals fj!04|) 
implies that q% = O^^). Therefore u can be neglected wrt 77, and the 
coefficients q\ , n > 2 can be neglected with respect to u. 
Consequently 

u ~ h k f(W^(fi) + aW£>{fi))dn ~ i / R /(-ra/x 2 + £ - arf - 2 M r 3 ) 

(HI) 

and formula (1 107j) is derived. Known the first breaking time from the well 
known formula 

n = TF ^ 7 = mm cmTF ^, F'(( b )<0, (112) 

equation (11061) implies that, if a > 0, the first breaking takes place when 
tf, = — 00 at yi, = 00. As we shall see now, to have breaking in a finite point 
of the (x, y, t) space, we need to consider n even. 

If n = 4, then 

W(n) = -tfi 2 + £-2ut + a (V^ 4 ) = 

-t/i 2 + +a[// - 4r?/i 3 ] - (r 4 - Au)/j, 2 + ( 113 ) 
a(-4r/ 3 - 8r?w + 4^ 2) )/i + £ - 2wr 4 + a[(r/ 4 - 8^gf ) + 4^ 3) + 6m 2 )], 

where 

r 4 = t - Qarf. (114) 

Since now W depends on the three unknowns u, n — 2, 3, the system of 
three equations (11041) . for n = 1,2,3, is an algebraic system characterizing 
a family of implicit solutions of dKP parametrized by the arbitrary real 
function / on a single variable. 

Following the same derivation as in the previous example, one obtains 
the following simple one dimensional implicit asymptotic formula for such a 
family of solutions 

u~^F(t + ar]±-2un), 

77 = 0(1), i + ari A -2un = 0(1), * > 1, < a < 1, 1 ; 

where F is again defined by 
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Known the first breaking time r& from (jll2p . equation (11141) implies that, 
if a > 0, the first (physical) breaking time t b is achieved at = (r] b = 
Ub/^tb = 0) and coincides with Tb, while Xb follows from 

Xb = ( b + 2F((b)Vt~ b . (116) 

We remark that formula fl 1 1 5 [) is an interesting particular case of the generic 
asymptotic formula (|35|) obtained in [33]- Summarizing, the first wave break- 
ing takes place at tb = u, where Tb and (b are given in (11121) . at the point 
{(b + 2F(Cb)Vh, 0) of the {x, y) plane. 

Now we are ready to write the case of an arbitrary power n. In this case 

W(fi) = -t/j, 2 + £-2ut + a (V^ n ) . (117) 

Since now W depends on the (n — 1) unknowns w, n = 2, . . . , n — 1, 
the system of {n — 1) equations (11041) . for n = 1, . . . , n — 1, is an algebraic 
system characterizing a family of implicit solutions of dKP parametrized by 
the arbitrary real function / on a single variable. 

Proceeding as in the previous examples, one obtains the following simple 
one dimensional implicit asymptotic formula for such a family of solutions 

M ~ i F(£ + (-l) n arf -2uT n ), , . 

V = 0{1), £+ {-l) n a V n -2uT n = 0(1), t > 1, < a < 1, v ; 

where 

rn = t-(-l)"a( n ^ 2 )^- 2 (119) 

and F is defined by fl89l) . Following the same considerations made in the 
particular case n = 4, if n is even and a > 0, the localized solution breaks 
at time t b = T b , where Tb is defined in (I112p . in the point (a^O), where 
Xb is defined in (11161) . We remark that also formula f 1 1 1 8 j) is an interesting 
particular case of the generic asymptotic formula (135|) obtained in [33]. 

4 Solvable NRH problems for the heavenly 
equation 

In this section we consider two examples of solvable NRH problems allow- 
ing to construct two different classes of similarity solutions of the heavenly 
equation. 
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4.1 Example 1. The invariant ipi^ an( ^ similarity so- 
lutions |47] 

If S(q,p) = qp, then, as in §3.1, the NRH problem reads 

^+ = ^ e i W^^) i ^+ = ^- e -*/Wr^. A ) (120) 

satisfying the symplectic and reality constraints, where now / is an arbitrary 
real function of two variables. Then the invariance equation (j52p becomes: 

= fo^- = zt\ 2 - (xt + yz)X + xy- z6 x + t6 y = W(\), (121) 

and the NRH problem linearizes and decouples: 

= ^- e V(w(A),A) ) ^+ = ^- e -if(W(^), (122) 

Equations (155)) become 

^+ e i(-y/+W = ^-eH-Vf-W = Ej(X), j = 1, 2, (123) 

where 

/ ± W S S/.FM/TO.A'). (124) 

It follows that 

(A) = x-zX-z<f> 
E 2 (X) = y-tX + t< f >, 

implying the following explicit solution of the NRH problem: 



(125) 



i(>f = E j (\)e i (-> i+1 f + M, j = 1,2. (126) 

A characterization of the corresponding solutions of dKP is obtained ob- 
serving that, since W in ( 11211) depends on the unknowns 6 x ,6 y , isolating 
the 1/A terms of equations fl 1 2 6 [) for |A| 3> 1 and comparing them with the 
asymptotics (1151) . we obtain the algebraic system 

9 y = x < f > -z(< A/>+f</> 2 ), 

^ = Z/</>-t(<A/>-i</> 2 ) ^ 

for the unknowns 6 x ,6 y . The constructed solutions of heavenly, parametrized 
by an arbitrary real function of two arguments, correspond to the following 
differential reduction: 

- W-! = t0 t + y9 y - x6 x - z9 z = 0. (128) 
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Substituting its general solution 

6 = A(x, y, z), x = xt, y — y/t, z — zt, (129) 

into the heavenly equation, one obtains the following similarity reduction of 
O): 

xA xx — yA%y + zAxz — Ayz + AxxAyy — A xy - = 0. (130) 

Therefore the algebraic system fl 1 2 7[) characterizes the above similarity solu- 
tions of heavenly. 

4.2 Example 2. The invariant (i/^i) 2 + (V^) 2 an d rota- 
tionally invariant solutions 

If S(q,p) = (q 2 +p 2 )/2, then equation (HTj) reads 



q(t) = cosH'(S)(t - t ) q + smH'(S)(t - t ) p , 
pit) = - sin H'{S) {t - t ) q + cos H\S) (t - t ) p , 

corresponding to the NRH problem 

ijf = cosh/ \) ^ + ismhf ( ^ 2+ 2 ^ 2 ,\) ^ 



(131) 



2 ) /y J Yl -Tboxxiuj I 2 J /x J Y2 ■ 



(132) 



-2 



satisfying the symplectic and reality constraints. Then the invariance equa- 
tion f )52|) becomes: 



-2 ,+2 ,_2 ,_2 



= = ^A 2 - (a* + aty)\ + ^+ (133) 

z9 y -at9 x = W(A), 



and the NRH problem linearizes: 



^ + = cosh(/(H/,A)) Vi _ -isinh(/(W,A)) ^ 2 ", 
^+ = isinh(/(W,A)) VT + cosh(/(W,A)) 



Equations become 



cosh /^(Xjipf - i sinh ^(X)^ = E^X), 
i sinh f ± (X)ipf + cosh / ± ( A)^ = £ 2 (A), 



(134) 



(135) 
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where 

E 1 (X) = -z\ + x-t(f), 
E 2 (X) = -t\ + y + z(f). 

Therefore the explicit solution of the NRH problem reads 



(136) 



^f(A) =cosh/ ± (A)£ 1 (A)+*sinh/ ± (A)£ 2 (A), 

^±(A) = -ismhf ± (X)E 1 (X) + cosh / ± (A)E 2 (A), K } 

A characterization of the corresponding solutions of heavenly is obtained 
observing that, since W in (11211) depends on the unknowns 6 x ,6 y , isolating 
the 1/A terms of equations f)126p for |A| ^> 1 and comparing them with the 



asymptotics (1431) . we obtain the algebraic system 

6 y =y<f>-t(<\f>-l <f> 2 ), , . 

9 X = x < f > -z (< \f > -§ < f > 2 ) [L6H) 

for the unknowns 9 X , 6 y . These solutions, parametrized by an arbitrary real 
function of two arguments, correspond to the following differential reduction: 

- = z6 t + xO y - ay6 x - atO z = 0. (139) 

Substituting its general solution 

9 = B(a, a = x 2 + ay 2 , (3 = z 2 + at 2 , 7 = xt — yz, (140) 

into the heavenly equation, one obtains the following reduction of (J5]): 



25 7 + 75 77 + 2aB a ^ + 2(3 B M + Aa^B a ^ + 2B a ((3B^ + 2aB a + 
AaaB aa + 4a 7j B a7 ) + 4(a/3 - ~? 2 )(B aa B„ n - B 2 ai ) = 0. 



(141) 



Therefore the algebraic system (11381) characterizes a class of rotationally 
invariant solutions of heavenly. 
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